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Abstract 

In this work, we are interested on the implementation of single-qubit gates on 
coupled Bose-Einstein condensates (BECs). The system, a feasible candidate 
for a qubit, consists on condensed atoms of different hyperfine levels coupled 
by a two-photon transition. It is well established that the dynamics of coupled 
BECs can be described by the two-mode Hamiltonian which takes into account 
the confinement potential of the trap and the effects of collisions associated 
with each condensate. Other effects, as collisions between atoms belonging to 
different BECs and detuning are included in this approach. We demonstrate 
how to implement two types of quantum logic gates: transfer-population gates 
(NOT, Y, and Hadamard), which require inversion of population between hy¬ 
perfine levels, and phase gates {Z, S, and T), which require self-trapping. We 
also discuss the experimental feasibility, checking the robustness of quantum 
gates against variations of physical parameters out of the ideal conditions for 
implementation of each quantum logic gate. 

Keywords: Bose-Einstein Condensates, Quantum Information, Quantum 
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Quantum information processing (QI^ is one of the remarkable research 
topics in Physics during the last decades [l], Q . Polinomial versus superpolino- 
mial times in factorization and faster search algorithms ii are two examples 
of the advantages of this kind of processing. The first requirement for the im- 
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plementation of QIP is the definition of a quantum of information or qubit Q , 
which can be done using a two-level physical system. In addition to associate 
the logic states 0 and 1 with one of two levels, the dynamics of the system should 
allow the preparation of superpositions of these states, which must be robust 
against decoherence. Several systems has been pointed out as good candidates 
for the implementation of a qubit. Those include photons, ions, charge and 
spin in semiconductor nanostructures, JosOThson junctions of superconductors, 
nuclear spins, and trapped neutral atoms [2|. 

Concerning trapped atomSjthe Bose-Einstein condensates (BECs) with al- 
kalii atoms, created in 1995 [a, [7|, has been pointed out as promising candidates 
for implementation of QIP proposals because the long decoherence time scale 
and the high degree of manipulation and control of its physical properties. In 
particular, coupled BECs can be implemented considering two experimental se¬ 
tups. In the first one, atoms condensed in two different sites of an optical lattice 
are coupled by tunneling d, . The second configuration consists of condensed 
atoms in different hyperfine levels coupled by laser BEl- The two-mode 
model provides an adequate theoretical description of the dynamics of both sys¬ 
tems 12| . Experimental achievements on coupled BECs include the observation 
of coherent oscillations 
a chip 


0,i imled condensates inside 

3 3. the creation of atomic interferometers [3 j and entangled quan¬ 
tum states [3l- Because the process of transfer of population are performed 
with high quality in a full controlled experimental setup, the system can be 
considered, in first approximation, in an equilibrium state. Although, the effect 
of the presence of non-condensed atoms over the coupled condensate are one 
of the causes behind of the destruction of condensate, the time scales are long 
enough for QIP operations. 

There are few theoretical proposals for the implementation of gates in cou¬ 
pled BECs. Those include the works of Calarco et al. [313 , focussed on 
the analysis of the conditions for the implementation of a phase gate by ma¬ 
nipulation of confinement potentials and a quantum logic gates using marker 
atoms [iOl' Pachos and Knight 211 proposed two-qubits operations including a 
Toffoli gate by combining both, tunneling and adiabatic passage of condensate 
atoms in a superlattice. Lee et al. explore the dipole-dipole interaction as a 
coupling mechanism [31 ■ In a recent works, Byrnes et al. 3 241 propose the 
codification of a quantum bit using collective states of coupled BECs by using 
an atomic coherent state, as well as the implementation of the two-qubit gate 
CNOT and an application to Grover algorithm. The authors point out that this 
proposal brings a natural advantage once the energy scale of the interaction is 
boosted by a factor corresponding to the number of condensed atoms, with the 
consequence of reducing the time of CNOT gate by the same factor. Concerning 
single-qubit operations, it is discussed how to rotate the qubit around z. 

Nevertheless, the practical implementation of single-qubit quantum logic 
gates was not sufficiently explored. To fill up this lack of information, we present 
our detailed analysis of the dynamics of coupled BECs which demonstrates how 
to implement the whole set of single-qubit quantum logic gates. Our physical 
system consists of atoms in a magneto-optical trap, condensed in two different 
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hyperfine levels which are coupled by a two-photon transition. First, we obtain 
an analytical solution for the time-dependent Schrodinger equation, written as 
an evolution operator. Then, we demonstrate that dynamics can be rewritten 
as successive applications of rotation operators, Ri. Finally, we explore this 
general result in order to found the conditions over the physical parameters 
which control the dynamics of the system in such a way that corresponds to 
the action of six quantum logic operations: the gates NOT, Y, and Hadamard, 
which require transfer of population, and the phase gates Z (tt), S (7r/2), and 
T (7r/4), which requires the inhibition of transfer of population using strong 
detuning. 

Our system consists of condensed atoms in two different hyperfine levels 
labeled by a and b. The interactions between the atoms are well described by 
assuming two-body collisions and the coupling of the two hyperfine levels is 
done via a two-photon transition. The second quantized Hamiltonian is written 

as {h = 1 ) [n, [ii, m, [23 


H = UJafla + Cbbflb + Janl -f yfeUfo -I- 2jabnanb 




( 1 ) 


where fia = a, fib = 6m, and we define the following physical parameters 
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Here, uij describes the effect of harmonic potential over the condensed atoms, 
depends on the scattering length Aj , describing collisions between the atoms 
condensed in the same hyperfine level a or b, and jab is associated with the 
collisions between atoms in different levels. The coupling parameter g takes 
into account the effect of a two-photon coupling between two different hyperfine 
levels with H being the coupling strength between the modes. This coupling 
could be resonant (A = 0) or not (A ^ 0). In equation ([T]), we also use the 
auxiliary parameters uj(^a,b) = ^{a,b) - l(a,b)- 

The task here is to find the evolution operator P{t) which evolves some 
initial state. The evolved state, which is the solution of the Schrodinger equation 
associated with Hamiltonian o, can be written as 

\^{t)) = p{t)\^m. ( 3 ) 

To find P{t), we use a technique which takes advantage of the properties of 
unitary transformations. We apply successively two transformations, being the 
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first one, 


U{t) = 

which removes the explicit temporal dependence from original Hamiltonian O- 
The transformed Hamiltonian is non-diagonal in the Fock basis and nonlinear, 
once there is a term which depends on An^ operator. The coupling is given by 
the nonlinear parameter A defined as. 


A = i (7a + 7b - 27ab) • 


(4) 


In the case of Rubidium collisions it applies that A = 0 [ifl, [ll|, although 


the manipulation of this parameter is possible by using Feshbach resonances 
which affects collision between atoms from different hyperfine levels [29l| . 
Concerning the non-diagonal feature, a second transformation 

V = 

is used to obtain the final form of a transformed Hamiltonian H'^ = H^V 
given by 


= UJoN + -fabN^ + 


-|- u; 2 N^ cos ^ -I- g sin ^ 


An. 


(5) 


Here we define uiq = {uja + <^b) /2, wi = (wa — Qb — A) /2, ui 2 = (7a — 7b) /2 and 


^ = arctan 


( _ 2 _^ 

bJi+UJ2N J ■ 


N = ha + n-b is the number operator, which accounts the 


L01-\-UJ2N 

total number of bosons in the system. Still, because is a conserved quantity, 
it will works as a c-number. After all this considerations, we obtain the final 
form for the evolved state (jS]) which reads as 

1^- {t)) = U (t) |«'(0)). 

That means that P{t) is given by 


P{t) = U (t) He 


-iH'' t 


Ft. 


( 6 ) 


This evolution operator can be seen as a product of rotation operators, because 
of the connection between two quantum harmonic oscillators and pseudo-spin 
operators [s^, so Jx = a^h + ah\ Jy = —i — a&t^ , and Jz = a)a — b'^b. At 
the same time, the operators inside the general expression above are, in fact, 
rotation operators BUI, once the operators p (i = x, y, z) are generators of 
rotations around ith axis. The final form for P{t) as a function of rotation 
operators Ri is 


P (r?, A, ro, tt) = e ^'^*Rz (At) Ry {-() Rz (rot) Ry (0 
where the auxiliary parameters are defined as 

. f 25 

^ = arctan 


(7) 


77 = 


r- A, 

(r — A) cos ^ -I- 2g sin 

Wb) -I- (7a + 7b) {N - 1)] N. 




( 8 ) 
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Here we define the frequency-scattering detuning as 


r = Waf, + (7a - 7h) - 1) , 


( 9 ) 


being uiab = oJa — the detuning between the effective frequency trap for each 
condensate. The parameter F quantifies differences between the trap frequencies 
and collision parameters for each condensate. This quantity and the non-linear 
parameter, A, play important roles in the definition of the necessary conditions 
for the implementation of quantum logic gates. 

After the obtention of the evolution operator, we proceed to exploit the pos¬ 
sibilities of applications of coupled BECs within the context of quantum infor¬ 
mation processing. In the standard quantum computation approach, encoding 
information requires the preparation of an initial state, which is defined in the 
computational basis {|0),|1)} as a coherent superposition with coefficients a 
and /3. Definition of qubits using coupled BECs are found in the literature, 
particularly in the context of condensed atoms in a two-sites potentia l 13 ill . Our 
proposal is to encode the qubit using an atomic coherent state (ACS) [3^, which 
can be created with hyperfine coupled BECs in two steps: i) the first one is the 
application of a pulse with a given duration to transfer the population from 
one condensate to another. For instance, if a 7r/2 pulse is applied over a BEC 
containing all condensed atoms, half of them will be in each hyperfine BEC 
after the pulse; and ii) the relative phase between the coupled BECs is created 
through the free evolution of the system ll|, [ig, jSSj . The ACS is defined as 
follows 


= 


1 
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( 10 ) 


which can be written as a quantum bit defined in the Bloch sphere as, 

\e,(j)) =cos^|0)-hsin^e*'^|l) = a |0)-h/3 |1), (11) 

with |0) = |0, A^) and |1) = |A^, 0), 6 and (p being the polar and azimuthal 
angles on Bloch sphere, respectively. Here, the polar angle 9 is related with the 
difference between atomic population of hyperfine levels while (p is the relative 
phase between them. The state (Hoj describes a particle with spin N/2, whose 
normalized mean values of the operators are given by < Jx > /N = sin 0 cos 
< Jy > /N = sin 0 sine/), and < Jz > /N = cos 6. ACS is analogous to the 
coherent state of the harmonic oscillator: it is the result of the application 
of a displacement operation over the ground state \N/2,—N/2) = |1) and is 
a minimum uncertainty state which fulfills the Heisenberg angular-momentum 
uncertainty relation | ^ where the mean values 

are calculated in a rotated coordinates system so z' is an axis which pass through 
the center of the ACS [s^. 


5 











Because P{t) is a succession of rotations, an initial state |'I'(0)) prepared as 
an ACS evolves to another ACS with new coefficients a{t) and /3(t). To fulfill our 
goal of finding the experimental conditions for implementation of single-qubit 
quantun^ates, we use the 2x2 matrix representation of rotation operators Rz 
and Ry [l| and obtain the matrix form of P given by 


where 




Pii 

Pl 2 

P21 

P22 


( cos^ ^ sin^ i 


1 
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= 2i cos - sin - sin 
2 2 
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e 2 


n i i ■ f 

2i cos - sm - sm — e 2 

2 2 V 2 / 


e* e * ( sin^ | cos^ 


( 12 ) 


(13) 


Now we are able to find the necessary conditions for the realization of a spe¬ 
cific single-qubit gate. This operations can be sorted in two sets, depending on 
the type of dynamics required for each gate. The first set, the transfer-population 
gates, contains the NOT, Y, and Hadamard gates and requires transfer of pop¬ 
ulation as well as changes on the relative phase. From the geometrical point 
of view, the evolution operator P{t) will change the polar angle and azimuthal 
angle of the qubit on Bloch sphere. The second set, the phase gates, given by 
the Z, S, and T gates, requires the inhibition of transfer of population, together 
with the imprinting of the required relative phase 0. 

To obtain the conditions for physical implementation, we proceed to com¬ 
pare the elements of evolution operator P, given by equation (1131) . with the 
corresponding elements on the matrix representation for each gate. Let us con¬ 
sider the NOT quantum gate. Comparing both, the matrix form for this gate 
and equation (ca, the elements of P must be given by Pn = P 22 = 0 and 
P 12 = P 21 = 1. After some algebra, we found that the detuning must follow 
the rule A = ^. In experiments, it is possible to control the detuning between 


^It is useful to remember the matrix form of gates NOT, Y and H, which are written as 




X = 


while the phase gates {Gp) are given by the general expression 


/ 1 0 \ 

VO J ’ 


with the value of relative phase ip corresponding to Gp = Z {ip = tt), S {ip = '7r/2), and T 
{(p = 7^1 A) quantum gates, respectively. 
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the two-photon transition, A, so we fix the gate evolution time as tNOT = 
Also, the values of the auxiliary parameters defined in equation (|8]) must fulfill 
the following conditions 

. f ‘2g \ tt 

t = arctan -- = —, 

^ \T-AJ 2 ’ 

tJ7 = (r — A) cos^ + 2gsin^ = 

V=\{{^a+^^b) + {la+7b){N -1)]N = ( 14 ) 

The solution of this coupled equations define two physical requirements: first, 
the value of the two-photon transition. A, which must be set as A = Aq = dg. 
Second, the frequency-scattering detuning, F, given by equation ([9]), must follow 
the condition F = Tq = Ag. 


Table 1: Evolution times and necessary conditions for the implementation of the transfer- 
population gates given by NOT, Y, and H. Here g is the coupling parameter, Aq is the 
detuning of the two-photon coupling, and Tq is the frequency-scattering detuning, as defined 
in equation Ml 


Gate 

^gate 

Ag 

Tg 

NOT 

ZTT 

A 

45 

4ff 

Y 

TT 

A 

25 

2<? 

H 

277 

A 

'^3 

—2g + A 


Following a similar procedure, we found the values for fgate, Aq, and Fg for 
all transfer-population gates, summarized in tabled] We note that the detuning 
of the two-photon transition, Aq, depends on coupling parameter g with the rule 
Aq = n^g, with being a different factor for each gate. The operation time is 
also a function of the coupling parameter g, which is responsible for the transfer 
of population, and all gates require a non-resonant two-photon coupling as well 
as specific non-zero values of the Fq parameter. This is explained because this 
two parameters (Aq and Fq) control the gain of relative phase between | 0 ) and 
|1). For coupled BECs on a chip [l5|, it is possible to control Uab while — 7 b 
is a fixed value. On the other hand, in the experimental setup in Ref. [l 6 | the 
values of collision parameters 7 ^ can be manipulated via Feschbach resonances, 
leaving uJab fixed. 

Concerning the phase gates, once that the goal is to gain relative phase 
suppressing the transfer of population, we expect that the requirements for 
the implementation of this gate includes a detuned two-photon transition. For 
example, in order to obtain the gate Z, it is necessary that the elements of 
propagator P follow the conditions Pu = 1, P 22 = — 1 with P 12 = P 21 = 0. 
This requirements are fulfilled when the evolved time is given by tz = 7 r/ 2 A. 
After some algebra, we obtain the conditions for the auxiliary parameters which 
are 

, f 2 g 

C = arctan ( ^ 
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Table 2: Evolution times and necessary conditions for the implementation of the quantum 
phase gates (Z, S, and T). Here g is the coupling parameter, Aq is the detuning of the 
two-photon coupling, and Eq is the frequency-scattering detuning, as defined in equation (0. 


Gate 

^gate 

Ag 

Tg 

Z 

TT 

2A 

^ 3 

-2A 

S 

Stt 

2A 

>3g 

A 

3 

f 

TT 

2A 

»4g 

A 

2 


tu = (r — A) cos^ + 2 ( 7 sin^ = 3A, 

^ ^ + Wb) - (7 q + 7b)] N + = A. 


(15) 


There are two possible mathematical solutions for the coupled equations. One 
solution is given by the condition g = 0, which is not interesting because de¬ 
scribes a physical situation where the condensates are uncoupled. The second 
solution is obtained when A = Aq ^ 2g/3, which implies a strong detuning 
of the two-photon transition as expected. For S and T, we obtain a similar 
condition for Aq and, once a value of Aq is fixed, a careful setting of quantity 
r = Tq is required in order to obtain the correct phase gain. With this con¬ 
siderations, we obtain the conditions for implementation of all phase gates that 
summarize in table 

It is interesting to check the dynamics by calculating the mean values of 
the operators Ji {i = x, y, z) and follow their evolution over the Bloch sphere as 
shown in figure[TJ First, we simulate the dynamics for |4'(0)) = a |0)-|-/3 |1) with 
a = cos (ff), and /3 = sin (|^) under the conditions for the implementation of 


gates NOT, figure[T]]a), F, hguredDb), and the phase gate Z, figure HDc). For 
all cases, the initial state evolves to the target final state, in perfect agreement 
with the action of each quantum gate operation: l'I'(tNOT)) = a\l) + P |0), for 
the NOT operation shown in figureUJa), |'I'(tY)) = crjl) —/3|0), for Y operation 
in figure dDb) , and |4'(tz)) = a |0) — /3|1) for Z operation in figure [TJc). The 
action of the Hadarmard gate is shown in figure (Ud) considering a different 
initial condition being |4'(0)) = |0). As can be seen, the gate creates the super¬ 


position 


|o)-Hi) 

V2 


We perform other simulations, not shown here, corroborating 


that initial state given by |1) evolves to the superposition , and the same 

gate takes both superpositions to the original initial states. From this results, 
our expectation is that the guidance provided by the tables [U and [21 together 
with a careful manipulation of physical parameters on coupled BECs, become 
a practical path to follow to guarantee the implementation of all single-qubit 
quantum gates. 

Concerning the feasibility of our proposal, the first task is to compare the 
evolution times for the quantum gates with the time scale of the typical deco¬ 
herence processes for BECs. Let us consider the value of the coupling parameter 
for experimental setup of Ref. , to calculate the evolution time fgate for each 
gate. To estimate the value of coupling g we use the reported value of t ^/2 = 170 
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of Ref. 0 , associated with the application of a 7r/2 pulse, which means the 
two-photon coupling is g = 27r x — = 27r x 1.5kHz. Using the table [U we 

obtain that the evolution times associated with NOT, Y, and H gates are given 
by Inot = ty = 0.16 ms, and = 0.12 ms, respectively. These values are 
significantly lower than the decoherence time, with value around 3 s (^ . For 
the experimental setup of Ref. [l^ , the value of the coupling parameter can be 
switched between zero and g = 2 tt x 600 Hz. Considering the highest value of 
g, we found the time scale for the gates as tNOT = ty = 0.42 ms, and % = 0.29 
ms which is 10“^ times smaller than decoherence time. Because the phase gates 
require a detuned interaction (Aq ^ g), the time scale goes as tgate oc I/Aq. 

The second and third tasks consist on check the effects of nonlinearity 
(A ^ 0), and the sensibility of dynamics against changes on values of frequency- 
scatering detuning, T, and the two-photon detuning, Aq, as dehned in tables [T] 
and [21 Hao and Gu [3^ studied the dynamic associated with Hamiltonian p]). 
considering non-zero values for nonlinear condition. The nonlinearity, which de¬ 
pends on J^, would affect our qubit, once the ACS becomes a squeezed state 0 
as time increases. This coupling also induces self-trapping, which might reduce 
the fidelity of quantum gates NOT, Y and H. Changes on specific conditions for 
r and Aq, could affect the fidelity of the implementation of quantum gates, be¬ 
cause they are crucial for both, the necessary gain of relative phase and transfer 
of population. To quantify the robustness of the gates, we perform simulations 
with the goal of solve the Schrodinger equation to calculate the exact state 
l^(^gate))- We compare the match between the target gate state. Ilk gate), and 
I'k (igate)), calculating the fidelity defined as 

= K'l^gate I'k (tgate))| ■ (16) 

For our simulations we consider N = 1000, the initial state being \0,(j)) = 
Itt/S, Ob and values for physical parameters according with experimental setup of 
Ref. [l 6 | . In order to quantify the effects of nonlinearity, we run the numerical 
simulations increasing the value of 7 ab in order to obtain A 7 ^ 0. Using a 
similar procedure, we quantify the effects of F parameter defined for each gate 
on tables d] and 121 considering increasing values of ujab so AT = |Fs — FgI ^0, 
being F^ the value used at simulation and Fq the value in tables |T] and |21 The 
effect of changes on A condition is quantified by considering SA = jA^ — Ag|, 
being A^ the value used at simulation. 

In figure 121 we plot our results of fidelity as a function of both, A and AF/F, 
considering four gates (a) H, (b) X and Y, {c)S, and (d) Z. For transfer- 
population gates, it is observed that the exact conditions for A and Fg obtained 
by analytical calculations are not exclusive for the implementation of a gate 
with high fidelity: the dark region with F = 1 extends to several choices for 
A and AF/F, with a linear relation between the two parameters for Hadamard 
gate, figure |2Ka), and a less extended region for X and Y gates, figure |2Db). 
This two behaviors is understood checking the definition of F parameter, equa¬ 
tion once this quantity depends on the collision parameters 7 ^, for small 
values of A the effect of nonlinearity can be compensated by modifying the 
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frequency-scattering detuning, which provides an extended set of experimental 
values with high fidelity gates. In the case of phase gates, the effect of AF over 
the decay of fidelity is weaker than the effect of the nonlinearity parameter. 
This is explained by the fact that all phase gates do not require a specific value 
of two-photon detuning, being enough to set this parameter as significatively 
higher than coupling parameter g. Once that the condition for Fq depends 
directly of A, it means that the variations around Fq ideal value on table [2] do 
not affect significatively the fidelity of gates. 

The last issue is to quantify the robustness of transfer-population gates out 
of the specific condition for detuning Aq defined on table [T] In figure O we 
show our numerical results for the three transfer-population gates. From our 
results, we observe that fidelity for all gates decays drastically when SA in¬ 
creases. Nevertheless, all three cases has the same behavior for small changes of 
A so 5 A/Ag < 0.1. For this interval of values, the fidelity has high values with 
F > 0.8, as indicated by the black dash arrow in figure |31 Additionally, the 
Hadamard and the Y gates continue with equal behavior until SA/Aq = 0.17, 
with F > 0.6. After that value, the X gate fidelity decays slower than the other 
two gates. The last gates have the same behavior until 5A/Aq = 0.15 (grey 
dash arrow in figure ED, with Hadamard being the less robust gate against small 
changes of the detuning value. From our simulations, we can conclude that fluc¬ 
tuations below 10% can be considered small, because they are associated with 
high values of fidelity. Because A affects both, the gain of relative phase and 
transfer of population, the effect of imprecision on the exact value of two-photon 
detuning at Aq increases when the gain of relative phase is mandatory on the 
implementation of the gate, which is the case for Y and H. 

Summarizing, we explore the properties of coupled Bose-Einstein conden¬ 
sates with the specific goal of implementation of controlled single-qubit quan¬ 
tum gates. We solve the Schrodinger equation to explore the dynamics in order 
to obtain the necessary physical conditions for the implementation of two sets 
of gates: those which require transfer of population (NOT, Y and ff), and the 
phase gates. We check the behavior of a Bloch vector defined using the mean 
values of pseudo-spin operators Ji with i = x,y,z setting the necessary condi¬ 
tions for each quantum logic gate. We conclude that the careful experimental 
control of the different detunings of the system (A, uJa — uib and ja — 7h) leads 
to the implementation of high fidelity quantum gates. The feasibility of our 
proposal was explored through numerical simulations, considering the effects 
of nonlinearity and deviations from the required value for frequency-scattering 
condition, Fq, as well as the two-photon detuning Aq. Our results show that it 
is possible to implement gates with high hdelity out of the ideal conditions. An 
open question is the scalability of this proposal. One option is the use of optical 
lattices H, H, 4^, ^|, with each site containing this two-modes coupling 
Bose-Einstein condensates. An appealing path is to explore the possibility of 
the use of the quantum bus in Ref. [2^ , an extra coupling between two of these 
coupled BECs using a cavity. For both alternatives, it is necessary to evalu¬ 
ate how the extra coupling affects the gates which demands heavy numerical 
simulations. This problem are topic for future research. 
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Figure 1: Dynamics over the Bloch state considering the initial state |’F(0)) = q;| 0) + /3|1) 
under the action of (a) NOT, (b) Y and (c) Z considering a = cos and /3 = sin (ff-)- 

(d) Action of the Hadamard gate over the initial state |’F(0)) = |0) The blue arrow denote 

the initial state, while the red one denote the final state and (x^y^z) —>■ ^ ^ 
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Figure 2: (Color online) Fidelity of {a)H, (b) X and Y, (c) S, and (d)Z gates as a function 
of A and AF/F 
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Figure 3: (Color online) Fidelity behavior for X (black solid line), Y (red triangles) and H 
(blue open squares) gates as a function of Aq- 
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